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ABSTRACT 

Aims. To determine the solar transition region and coronal radius at EUV wavelengths and its time evolution during 
Solar Cycle XXIII. 

Methods. We use daily 30.4 and 17.1 nm images obtained by the Extreme Ultraviolet Imager (EIT) aboard the SoHO 
satellite and derive the solar radius by fitting a circle to the limb brightness ring. 

Results. The weighted mean of the temporal series gives (967". 56 ± 0".04) and (969". 54 ± 0".02) at 30.4 and 17.1 nm 
respectively. No significant correlation was found with the solar cycle at any of the two wavelengths. 
Conclusions. Since the temperature formation of the 30.4 nm line is between (60 — 80) 10 3 K (Transition Region), the 
obtained result is bigger than that derived from present optical atmospheric models. On the contrary this height is 
compatible with radio models. 
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1. Introduction 

The solar atmospheric levels at which most of the different 
wave-band brightness are produced have been analyzed 
' by many authors to understand how the energy flow 
may affect the local equilibrium. Although there is a 
global hydrostatic equilibrium produced by an almost 
stable photospheric level, t he total s olar irradiance (the 
solar constant) varies (e.g. ISthd [20021 ). The constancy or 
variation of the photospheric solar diameter has been of 
a much debate subject. The coronal heating, however, 
is an indication that some energy is being deposited 
above the photosphere modifying the equilibrium. The 
chromosphere and transition region height variations 
have been verified by solar radius measurements at radio 
waveleng ths with clear cor relation with the solar constant 
changes (jCosta et alJll999f >. 

In recent years lEmilio et al.l (|2000l ) and iKuhn et al.l 
(2004J) have shown that the photospheric radius variation, 
if it exist, must be smaller than 15 mas during a whole 
solar cycle, therefore it can be taken as constant. Analysis 
in other frequencies of the electromagnetic spectrum 
have been carried out mainly in the radio domain. 
ISelhorst et al.1 (|2004l ) have reviewed previous results and 
determined the solar radius variation at 17 GHz using 
maps from the Nobeyama Radioheliograph (NoRH) with 
a space resolution of 5". They found a good correlation 
between the mean radius and the sunspot number but an 
anti-correlation between the polar radius and the sunspot 
number. 

At UV frequencies there are very few previous works, 
with short temporal ser i es an alysis. In one of these very 
few works, IZhang et alj (|1998[ ) used EIT full bun images 
and concluded that: i) the solar disk is prolate (the same 



conclusion was obtained by lAuchere et all Il998f) and ii) 
the chromospheric 30.4 nm limb is significantly higher 
than the limb at 17.1, 19.5 and 28.4 nm wavelengths, 
arguing that spicules can contribute to the increase of the 
chromospheric / Transition Region (TR) height. 

The determination of the solar atmosphere layer's 
height is relevant to the formulation of atmospheric 
models. Many of these models are based on optical spectral 
analysis. The density and temperature dependence on 
height are derived b y fitting a pred icted spect rum to 
the observed one, e.g. [Vernazza et ail (119731. Il976t Il98ll ); 
iFontenla et all (|1990l 119911 . 119931 l2002h . All of these mod- 
els obtain a TR height between 1700 and 2300 km above 
the photosphere. On the contrary, some Ha observations 
show a TR at around 5000 km above Photosphere (e.g. 
IZirinl Il99rfl . At microwaves, the solar radius has even 
bigger values ~ 10 4 km (e.g. I Selhorst et al .1120041 ). IZirinl 
(|l996l ) argued that the hydrostatic assumption of the 
atmospheric models is responsible for the disagreement 
betwe en predictions and observations. On the other 
hand, ISelhorst et al.l (|2005l ) claim that spicules are suf- 
ficient to increase the height of the solar limb when a 
model with a brightness temperature T& ~ 11 000 K at 
17 GHz and plasma densities few times 10 10 cm -3 are used. 

In this work we report the determination of the chro- 
mospheric/TR and coronal solar radius and their varia- 
tion during the solar cycle XXIII using full Sun images 
of the Extreme Ultraviolet Imager Telescope (EIT) aboard 
the Solar Heliospheric Observatory (SoHO) obtained with 
the passband filters centered at 30.4 nm and 17.1 nm. 
Our determination is based on the limb bright ring seen 
in both wavelengths. Since the 30.4 nm He line forms at 
around few times 10 4 K we com pare our observations with 
the microwave determinations of lSelhorst et al. (2004). We 
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present the method and results in Sect.[2J while in Sect. [3] 
we analyze the obtained time series and discuss our results. 



2. The determination of the limb bright ring 

The EIT is one of th e 12 i nstruments aboard the SoHO 
satellite (iMoses et al.l Il997|) . It obtains images in four 
wavelengths by means of passband filters centered on 
the coronal lines A = 17.1 nm (Fe ix,x, T ~ 10 6 , K), 
A =19.5 nm (Fe XII, T ~ 1.4 10 6 K) and A = 28.4 nm 
(Fe xv, T ~ 2.1 10 6 K). The fourth wavelength range 
corresponds to a chromospheric / TR line, A = 30.4 nm 
(He II Ly a, T ~ 6-8 10 4 K). 

We choose images centered at the wavelengths A = 17.1 
nm and A = 30.4 nm. The first line corresponds to the 
blending of Fe X, A = 17. 106 nm and Fe ix, A = 17.056 nm 
(|Cowan fe Peacockl I1965D . Both lines are optically thin, 
therefore t hey are most inte nse where the column depth is 
bigger (e.g. IWithbroelll970[ h The maximum of the Fe ix,x 
emission is a pure coronal feature which can be used to 
determine a solar coronal radius that can be compared 
with radio continuum observations. Furthermore, since 
the emission is optically thin, changes in the coronal 
energy input may affect this radius and therefore this 
definition allows us to investigate the coronal heating 
during a solar cycle. On the other hand, the maximum 
of the radial profile is very sharp (see Figure \2§ making 
the radius determination very precise. Our criterion for 
this l ine is different to the one used by IZhang et al.l 
(1998) who determines the occulting lim b of the EUV 
iron l ines in a way which is equivalent to lAuchere et alj 
(|1998f ) limb definition as the inflection point of the 
radial profile. We note although, that this occulting or 
blocking limb corresponds to a Chromospheric layer which 
is investigated in this paper by means of the helium images. 

The formation of the He II line is more controversial. 
There is an ambiguity between the predicted and observed 
intensities and a discussion about its mechanism (e.g. 
Uordanl 119751 : lAndretta etaU l2003f) . In general the line 
is considered optically t hick and formed in a thin layer 
(jPietarila fe Judge] l2004h resulting also in a limb bright 
ring. Although this ring is not so intense and clearly de- 
fined as the coronal one, it can be seen in the EIT images 
during the whole solar cycle. The line is contaminated by 
the coronal A = 30.3 nm Si xi line, which may have a 
relative intensit y of up to 20% above the disk brightening 
(|Aucherd l2000f ). Since Si xi is a coronal line, we expect 
that it will increase the radius determination, although 
there is no easy way to quantify the effect. We believe that 
since it is a weak line it will produce a small shift masked 
by the rest of the uncertainties. 

We define the solar radius as the semi-distance between 
the diametrically opposed maxima of the limb bright points 
of the rings. This definition allows us to determine the 
radius with an automatic methodology which is explained 
below. 



2.1. Limb Extraction and Radius Determination 

Every image was corrected using standard procedures in- 
cluded in the Sola r Softw a re pa ckage for the EIT instru- 
ment. We use the ICannvl (|l986f ) algorithm to get a first 
approach to the limb ring, a circular fitting to this ring 
gives us a disk center determination better than the one 
found in the image header. Once the center is well deter- 
mined, radial profiles are extracted every half a degree in 
polar angle. The limb brightening of each profile is fitted 
to a special function f(r) that depends on the wavelength 
band: 

17.1 nm : f(r) = a D atan((ai — r) a^) + (13 + 

a 4 r aB exp(-|r - r m | a 6 ) 
30.4 nm : f(r) — a exp (— (r — r m ) 2 /2a\) + 02 + a 3 r + 

CI4 r 2 

Fe ix, x limb brightening is better represented by a cusp 
function while the arctangent is used to represent the 
transition between the scattered coronal light and the quiet 
sun. He 11 limb profiles do not have a clear peak and thus 
are not well represented by a cusp function, therefore we 
use a Gaussian to determine the mean brightest position of 
the limb. In this case a second degree polynomial represents 
the transition from the scattered light and the quiet sun. 
In both functions, and r m are the fitting parameters to 
be determined, the latter is the position of the limb max- 
imum. The fitting is carried out by means of the amoeba 
algorithm (|Press et al.l I2002T ). setting the optical radius 
as the lower limit of the fitting interval (see Figure for 
examples of the fittings). We iteratively fit a circle to the 
set of (x, y) ring positions discarding at every step those 
that are outside the band (R — 1.8 a, R + 1.8 a), with R the 
circle radius and a the standard deviation of the fitting. 
The iteration goes until it converges with all pairs inside 
the fitting band. To convert from the CCD pixel units to 
arc seconds ("), we use the pl ate scale determined durin g 
a 1999 Mercury transit by lAuchere fe Artznerl (|2004l ), 
namely 2.627 ± 0.001 arc sec pixel -1 . This figure was later 
confirmed with a new transit of Mercury in 2003 (Artzner, 
personal communication) . The resulting size is corrected to 
1 AU using the satellite distance found in the image header. 



2.2. Results 

Figure shows the time evolution of the solar 17.1 and 
30.4 nm radii. Heights arc relative to the photosphcric 
radius at 1 AU, R Q = 959".68 = 695489.2 km. Points 
represent individual measurements, continuous curves are 
28 days running means and the horizontal thick lines are 
the weighted means. There are some gaps in the data 
when the instrument stoped its operation. Remarkable 
oscillations are observed before the gap at the end of 1998. 
An incorrect determination of the satellite distance, or 
a variation of the telescope focal distance produced by 
thermal oscillations can account for an incomplete data 
correction. Nonetheless it is not clear why the effect is 
so important at 30.4 nm and less significant at 17.1 nm. 
Since we do not have certainty on the origin of the strong 
oscillations we decide not to use the data for the analysis 
before the long observing stop of the instrument in 1998. 
In general, the mean uncertainty in radius determination 
is 1".7, almost half a pixel, at 30.4 nm and 1", a third of 
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Fig. 1. EUV limb brightening. Images show the limb during different solar cycle phases. Right panels represent the mean 
limb profile. 
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Fig. 2. Examples of limb fitting. Left, a limb profile at 17.1 
nm (gray) and the fitting (black). Right, the same for the 
30.4 nm limb profile. 
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Fig. 3. Time evolution of the solar radii. Points represent 
individual measurements. Solids curves are 28 days running 
means at 17.1 nm (upper) and 30.4 nm (lower). Straight 
thick lines are the weighted mean values for the period 1999 
- 2005. The vertical bars represent the mean uncertainty for 
one measurement. 



a pixel, at 17.1 nm. The bigger uncertainty for the line of 
He II reveals its weaker intensity and its complex nature. 

The very interesting harmonic variations of the solar 
limb during the cycle shown in [3] are still under analyses to 
be published somewhere else. 



3. Analysis and Discussion 

Table [1] summarizes our results. First we note that the 
weighted means show a clear different height for the chro- 
mospheric He II and the coronal Fe ix,x lines. Even if we 
use the mean uncertainty of one measurement the heights 
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Radius 

(") 



Height 

km 



He n 


967.56 


± 


0.04 


5710 


±30 


Fe IX, X 


969.54 


± 


0.02 


7146 


± 15 



Table 1. Solar radius and height over the photosphere. 



do not overlap. It is also clear that the chromosphere 
height is well above where present atmospheric models 
predict and near where iJohannesson fe Zirinl ( 19961) found 
it usin g Ha filtergrams and is comparable to IZhang et al.l 



lAuchere et al.l (|l998l ) have shown that the Sun is 
1".5 and 5" prola te at 17.1 and 30.4 nm respectively. 
Zha ng et al.l (|1998f) have also observed a similar prolat- 
ness of the EUV solar limb. Due to the noisy limb we 
did not fit an ellipse. We use instead a circle fit much 
simpler and faster. It should be noted that the ellipsoidal 
form is revealed aft er the average of over 400 images in 
lAuchere et al.l (|1998( ) work. Therefore, our results represent 
an averaged value over the solar limb. Taking Auchere's 
results into account we can correct our determination 
of the equatorial radius subtracting 2" and the polar 
radius adding 3". In that way we get 4130 km equa torial 
height and 7750 km polar height. I Zhang etTa l. ( 1998) find 
3100 ± 1200 km and 6600 ± 1200 km for the equatorial 
and polar heights respectively, in accordance with our 
results. On the other hand, we cannot compare our 17.1 
nm height with that obtained in Zhang's work because of 
our different definition of coronal radius. 



3.1. Comparison with radio observations and solar indices 
variations 



The S SC atmospheric model proposed by ISelhorst et al.l 
(|2005h predicts a solar radius of 966". 5 at 17 GHz. 
This microwave continuum emission comes from a solar 
atmospheric layer at a temperature of around 10 4 K, i.e. 
the chromosphere below the region where the Fe ix,x and 
He II are formed. Thus, the re sults of this wo r k agre e with 
that model. Not withstand, ISelhorst et alj (|2004l ) have 
determin ed a radius of aroun d 975" from 17 GHz NoRH 
maps. In S elhorst et al.l (|2005| ) the authors claim that the 
presence of spicules can explain the higher observed radius. 
Since typical s picules tempe rature is between 0.5 10 4 and 
2.5 10 4 K (see ISterlinel 120001 and references therein), their 
contribution to the EUV lines can be neglected and there 
is no contradiction in having the height of the 30.4 nm 
chromospheric line forming below the observed 17 GHz 
continuum. 

In Figure 2] we plot our results together with the 
daily Wolf Number and the daily mean flux density at 
10.7 cm solar indices. To ease the comparison we have 
rescaled the values in the band [—0.5,0.5]. The figure 
shows no clear indication of a close association between 
the radius and the solar indices. The cross correlation 
between the He n and Fe ix,x radius and the solar indices 
are shown in table [2j The near zero figures indicate that 
no correlation exists, although one is tempted to affirm 
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Fig. 4. EUV radius variations compared with solar cycle in- 
dices. From top to bottom: 17.1 nm radius, 30.4 nm radius, 
Wolf Number and flux density at 10.7 cm. All values are 
rescaled between [—0.5, 0.5] to facilitate the comparison. 





17 GHz 


Irr. 


Wolf N. 


10.7 cm 


He II 


+0.14 


+0.20 


+0.30 


-0.12 


Fe ix, X 


-0.12 


-0.38 


-0.43 


-0.40 



Table 2. Cross correlation between our results of the solar 
radii and solar activity indicators. From left to right: solar 
radius at 17 GHz, Total Solar Irradiance, Wolf Number and 
flux density daily means at 10.7 cm (S-componcnt). 



that there is some anti correlation, of around -0.4, between 
the coronal line and the solar indices. The real nature 
of such an anti correlation remains to be confirmed. In 
any case, it is surprising that the solar activity does not 
change the solar EUV radius neither at chromospheric 
nor at coronal heig hts as it is observed, v.g. at 17 GHz 
([Selhorst et alj 120051 ). As an explanation we think that 
the microwave continuum comes from a wide atmospheric 
layer receiving an important contribution from the spicules 
while a narrower layer is expected to form the EUV lines. 
Therefore it is possible that the variation in position 
of the EUV layers are not detectable with the present 
methodology and instrumentation. 

Although we find no long term variations and the use 
of the running means erases changes with durations of the 
order or lower than 28 days, Figure [3] shows fluctuations 
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with hundred days temporal scale, more sensitive for the 
He II than for the Fe ix,x images. This variation deserves a 
careful examination and a future study will have to show if 
it is just an instrumental effect or a real solar phenomenon. 
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